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ABSTRACT 

This paper makes the first systematic attempt to determine using perturbation theory 
the positions of images by gravitational lensing due to arbitrary number of coplanar 
masses without any symmetry on a plane, as a function of lens and source parameters. 
We present a method of Taylor-series expansion to solve the lens equation under a small 
mass-ratio approximation. First, we investigate perturbative structures of a single- 
complex-variable polynomial, which has been commonly used. Perturbative roots are 
found. Some roots represent positions of lensed images, while the others are unphysical 
because they do not satisfy the lens equation. This is consistent with a fact that the 
degree of the polynomial, namely the number of zeros, exceeds the maximum number 
of lensed images if N=3 (or more). The theorem never tells which roots are physical (or 
unphysical). In this paper, unphysical ones are identified. Secondly, to avoid unphysical 
roots, we re-examine the lens equation. The advantage of our method is that it allows 
a systematic iterative analysis. We determine image positions for binary lens systems 
up to the third order in mass ratios and for arbitrary N point masses up to the second 
order. This clarifies the dependence on parameters. Thirdly, the number of the images 
that admit a small mass-ratio limit is less than the maximum number. It is suggested 
that positions of extra images could not be expressed as Maclaurin series in mass 
ratios. Magnifications are finally discussed. 
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1 INTRODUCTION 

Gravitational lensing has become one of important subjects 
in modern astronomy and cosmology (e.g., Schneider 2006, 
Weinberg 2008). It has many applications as gravitational 
telescopes in various fields ranging from extra-solar planets 
to dark matter and dark energy at cosmological scales. This 
paper focuses on gravitational lensing due to a N-point mass 
system. Indeed it is a challenging problem to express the 
image positions as functions of lens and source parameters. 
There are several motivations. One is that gravitational lens- 
ing offers a tool of discoveries and measurements of plane- 
tary systems (Schneider and Weiss 1986, Mao and Paczynski 
1991, Gould and Loeb 1992, Bond et al. 2004, Beaulieu et al. 
2006), compact stars, or a cluster of dark objects, which are 
difficult to probe with other methods. Gaudi et al. (2008) 
have recently found an analogy of the Sun- Jupiter-Saturn 
system by lensing. Another motivation is theoretically ori- 
ented. One may be tempted to pursue a transit between a 
particle method and a fluid (mean fleld) one. For microlens- 



ing studies, particle methods are employed, because the sys- 
tems consist of stars, planets or MACHOs. In cosmological 
lensing, on the other hand, light propagation is considered 
for the gravitational field produced by inhomogeneities of 
cosmic fluids, say galaxies or large scale structures of the 
universe (e.g., Refregier 2003 for a review). It seems natu- 
ral, though no explicit proof has been given, that observed 
quantities computed by continuum fluid methods will agree 
with those by discrete particle ones in the limit N —> oo, at 
least on average, where A'' is the number of particles. 

Related with the problems mentioned above, we should 
note an astronomically important effect caused by the finite- 
ness of A'^. For most of cosmological gravitational lenses 
(both of strong and weak ones), a continuum approxima- 
tion can be safe and has worked well. There exists an excep- 
tional case, however, for which discreteness becomes impor- 
tant. One example is a quasar microlensing due to a point- 
like lens object, which is possibly a star in a host galaxy 
(for an extensive review, Wambsganss 2006). A galaxy con- 
sists of very large number particles, and light rays from 
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an object at cosmological distance may have a chance to 
pass very near one of the point masses. As a consequence of 
finite- AT effect in large N point lenses, anomalous changes in 
the light curve are observed. For such a quasar microlensing, 
hybrid approaches arc usually employed, where particles arc 
located in a smooth gravitational field representing a host 
galaxy. It is thus likely that N point-mass approach will be 
useful also when we study such a finite- iV effect at a certain 
transit stage between a particles system and a smooth one. 
Along this course, An (2007) investigated a point lens 
model, which represents a very special configuration that 
every point masses are located on regular grid points. 

For a A'' point-mass lens at a general configuration, very 
few things arc known in spite of many efforts. Among known 
ones is the maximum number of images lensed by A'^ point 
masses. After direct calculations by Witt (1990) and Mao, 
Fetters and Witt (1997), a careful study by Rhie (2001 for 
N=4, 2003 for general N) revealed that it is possible to ob- 
tain the maximum number of images as 5(A'' — 1). This the- 
orem for polynomials has been extended to a more general 
case including rational functions by Khavinson and Neu- 
mann (2006). (See Khavinson and Neumann 2008 for an el- 
egant review on a connection between the gravitational lens 
theory and the algebra, especially the fundamental theorem 
of algebra, and its extension to rational functions). 
Theorem (Khavinson and Neumann 2006): 
Let r{z) = p{z)/q{z), where p and q are relatively prime 
polynomials in z, and lot n bo the degree of r. If n > 1, 
then the number of zeros for r{z) — z* < 5(n — 1). Here, z 
and z* denote a complex number and its complex conjugate, 
respectively. 

Furthermore, Bayer, Dyer and Giang (2006) showed 
that in a configuration of point masses, replacing one of the 
point deflectors by a spherically symmetric distributed mass 
only introduces one extra image. Hence they found that the 
maximum number of images due to A'^ distributed lensing 
objects located on a plane is 6(A'^ — 1) -I- 1. 

Global properties such as lower bounds on the number 
of images are also discussed in Fetters, Levine and Wambs- 
ganss (2001) and references therein. 

In spite of many efforts on A'^ lensing objects, functions 
for image positions are still unknown even for N point-mass 
lenses in a general configuration under the thin lens approx- 
imation. Hence it is a challenging problem to express the 
image positions as functions of lens and source locations. 
Once such an expression is known, one can immediately ob- 
tain magnifications via computing the Jacobian of the lens 
mapping (Schneider et al. 1992). 

Only for a very few cases such as a single point mass 
and a singular isothermal ellipsoid, the lens equation can bo 
solved by hand and image positions are known, because the 
lens equation becomes a quadratic or fourth-order one (For 
a singular isothermal ellipsoid, Asada et al. 2003). For the 
binary lens system, the lens equation has the degree of five 
in a complex variable (Witt 1990). It has the same degree 
also in a real variable (Asada 2002a, Asada et al. 2004) . This 
improvement is not trivial because a complex variable brings 
two degrees of freedom. This siuglc-real-variable polynomial 
has advantages. For instance, the rmmber of real roots (with 
vanishing imaginary parts) corresponds to that of lensed im- 
ages. The analytic expression of the caustic, where the num- 
ber of images changes, is obtained by the fifth-order poly- 



nomial (Asada et al. 2002c). Galois showed, however, that 
the fifth-order and higher polynomials cannot be solved al- 
gebraically (van der Waerden 1966). Hence, no formula for 
the quintic equation is known. For this reason, some numer- 
ical implementation is required to find out image positions 
(and magnifications) for the binary gravitational lens for a 
general position of the source. Only for special cases of the 
source at a symmetric location such as on-axis sources, the 
lens equation can be solved by hand and image positions 
arc thus known (Schneider and Weiss 1986). For a weak 
field region, some perturbative solutions for the binary lens 
have been found (Bozza 1999, Asada 2002b), for instance 
in order to discuss astrometric lensing, which is caused by 
the image ccntroid shifts (for a single mass, Miyamoto and 
Yoshii 1995, Walker 1995; for a binary lens, Safizadeh et al. 
1999, Jeong et al. 1999, Asada 2002b). 

If the number of point masses A^ is larger than two, the 
basic equation is much more highly non-linear so that the 
lens equation can be solved only by numerical methods. As 
a result, observational properties such as magnifications and 
image separations have been investigated so far numerically 
for A'^ point-mass lenses. This makes it difficult to investigate 
the dependence of observational quantities on lens parame- 
ters. 

This paper is the first attempt to seek an analytic ex- 
pression of image positions without assuming any special 
symmetry. For this purpose, we shall present a method of 
Taylor-series expansion to solve the lens equation for N 
point-mass lens systems. Our method allows a systematic 
iterative analysis as shown later. 

Under three assumptions of weak gravitational fields, 
thin lenses and small deflection angles, gravitational lensing 
is usually described as a mapping from the lens plane onto 
the source plane (Schneider et al. 1992). Bourassa and Kan- 
towski (1973, 1975) introduced a complex notation to de- 
scribe gravitational lensing. Their notation was exclusively 
used to describe lenses with elliptical or spheroidal symme- 
try (Borgeest 1983, Bray 1984, Schramm 1990). For AT point 
lenses, Witt (1990) succeeded in recasting the lens equation 
into a single-complex-variable polynomial. This is in an el- 
egant form and thus has been often used in investigations 
of point-mass lenses. An advantage in the single-complex- 
variable formulation is that we can use some mathemati- 
cal tools applicable to complex-analytic functions, especially 
polynomials (Witt 1993, Witt and Fetters 1993, Witt and 
Mao 1995). One tool is the fundamental theorem of algebra: 
Every non-constant single- variable polynomial with complex 
coefficients has at least one complex root. This is also stated 
as: every non-zero single-variable polynomial, with complex 
coefficients, has exactly as many complex roots as its degree, 
if each root is counted as many times as its multiplicity. On 
the other hand, in the original form of the lens equation, one 
can hardly count up the number of images because of non- 
linearly coupled properties. This theorem, therefore, raises 
a problem in gravitational lensing. The single- variable poly- 
nomial due to A'^ point lenses has the degree of + 1, 
though the max;imum number of images is 5(A'^ — 1). This 
means that unphysical roots are included in the polynomial 
(for detailed discussions on the disappearance and appear- 
ance of images near fold and cusp caustics for general lens 
systems, see also Fetters, Levine and Wambsganss (2001) 
and references therein). First, we thus investigate explicitly 



Gravitational Lensing by N Point Mass 3 



behaviors of roots for the polynomial lens equation from 
the viewpoint of perturbations. We shall identify unphysical 
roots. Secondly, we shall re-examine the lens equation, so 
that the appearance of unphysical roots can be avoided. 

This paper is organised as follows. In Section 2, the 
complex description of gravitational lensing is briefly sum- 
marised. The lens equation is embedded into a single- 
complex-variable polynomial in Section 3. Perturbative 
roots for the complex polynomial are presented for binary 
and triple systems in sections 4 and 5, respectively. They are 
extended to a case of A'^ point lenses in section 6. In section 7, 
we re-examine the lens equation in a dual-complex- variables 
formalism and its perturbation scheme for a binary lens for 
its simplicity. The perturbation scheme is extended to a A'' 
point lens system in section 8. Section 9 is devoted to the 
conclusion. 



2 POLYNOMIAL FORMALISM USING 
COMPLEX VARIABLES 

We consider a lens system with N point masses. The mass 
and two-dimensional location of each body is denoted as Mi 
and the vector Ei, respectively. For the later convenience, 
let us define the angular size of the Einstein ring as 



(1) 



where G is the gravitational constant, c is the light speed, 

Mtot is the total mass Yl't-i ^^'^ ^L' ^S ^^'^ -^LS 
note distances between the observer and the lens, between 
the observer and the source, and between the lens and the 
source, respectively. In the unit normalised by the angular 
size of the Einstein ring, the lens equation becomes 



f3 = e 



e 



e 



(2) 



where (3 = {Px,l3y) and 9 — {9x,0y) denote the vectors for 
the position of the source and image, respectively and we 
defined the mass ratio and the angular separation vector as 
Vi = Mi/Mtot and = Ei/6E = {eix,eiy)- 

In a formalism based on complex variables, two- 
dimensional vectors for the source, lens and image positions 
are denoted as w = P^+iPy, z = 0^+i6y, and — Ei^+ieiy, 
respectively (See also Fig. [T]). By employing this formalism, 
the lens equation is rewritten as 



JV 



(3) 



where the asterisk * means the complex conjugate. The lens 
equation is non-analytic because it contains both z and z* . 



3 EMBEDDING THE LENS EQUATION INTO 
AN ANALYTIC POLYNOMIAL 

The complex conjugate of Eq. ([Sjl is expressed as 



Z — f 



(4) 




Obs 



Figure 1. Notation: The source and image positions on complex 
planes are denoted by w (the circle) and z (the filled disk), re- 
spectively. Locations of N point masses are denoted by (filled 
triangles) for i = 1, ■ ■ ■ , N . Here, we assume the thin lens approx- 
imation. 



This expression can be substituted into z* in Eq. ([Sjl to 
eliminate the complex variable z* . As a result, we obtain 
a (A*'^ + l)-th order analytic polynomial equation as (Witt 
1990) 



fc=i j^k 



e=i \ fc=i 

JV N 

JV / JV 



N N 



fc = l jT^fc 



(5) 



Equation (A3) in Witt (1990) takes a rather complicated 
form because of inclusion of nonzero shear 7 due to sur- 
rounding matter. Bayer et al. (2006) uses a complex formal- 
ism in order to discuss the maximum number of images in a 
configuration of point masses, by replacing one of point de- 
flectors by a spherically symmetric distributed mass. Their 
lens equation (3) agrees with Eq. ([Sj. In order to show this 



agreement, one may use ( — 1) 



JV+l 



It is worth- 



while to mention that Eq. ((5} contains not only all the solu- 
tions for the lens equation ([2]) but also unphysical false roots 
which do not satisfy Eq. ((2)1, in price of the manipulation 
for obtaining an analytic polynomial equation, as already 
pointed out by Rhie (2001, 2003) and Bayer et al. (2006). 
Such an inclusion of unphysical solutions can be easily un- 
derstood by remembering that we get unphysical roots as 
well as true ones if one takes a square of an equation includ- 
ing the square root. In fact, an analogous thing happens in 
another example of gravitational lenses such as an isother- 
mal ellipsoidal lens as a simple model of galaxies (Asada et 
al. 2003). 

In general, the mass ratio Vi satisfies < i/i < 1, so that 
it can be taken as an expansion parameter. Without loss of 
generality, we can assume that the first lens object is the 



4 H. Asada 



most massive, namely mi > rrii for i — 2,3, ■ ■ ■ , N . Thus, 
formal solutions are expressed in Taylor series as 



EE- 

P2=0P3=0 



E '^'^'^i 



(6) 



where the coefficients Z(p2)(p3) - (pjv) independent of Ui. 

Up to this point, the origin of the lens plane is arbitrary. 
In the following, the origin of the lens plane is chosen as the 
location of the mass mi, such that one can put ei = 0. 
This enables us to simplify some expressions and to easily 
understand their physical meanings, mostly because gravity 
is dominated by mi in most regions except for the vicinity 
of rrii {i 7^ 1). Namely, it is natural to treat our problem 
as perturbations around a single lens by mi (located at the 
origin of the coordinates). 

In numerical simulations or practical data analysis, 
however, one may use the coordinates in which the origin is 
not the location of mi . If one wishes to consider such a case 
of ei 7^ 0, one could make a translation by ei as z z + ei, 
10 — > ui + ei and + ei in our perturbative solutions 

that are given below. 



4 PERTURBATIVE SOLUTIONS FOR A 
POLYNOMIAL FORMALISM 1: BINARY 
LENS 

In this section, we investigate binary lenses explicitly up to 
the third order. This simple example may help us to un- 
derstand the structure of the perturbative solutions. For an 
arbitrary A'' case, expressions of iterative solutions are quite 
formal (See below). 

For simplicity, we denote our expansion parameter as 
m = 1/2- This means i^i = 1 — m. We also denote £2 simply 
by e. 

In powers of m, the polynomial equation is rewritten as 



^m''/fc(2) = 0, 

k = 



(7) 



where we defined 

ft)(z) = (z — e)^[(w* — e*)z + l\{w* z'^ — ww* z — w), 
fi{z) = (z-w) 

X (^e{z - w)[{2w* - €*)z + 2] - 6*2^(2 - e) - ezj , 

f2{z) = e^(z~w). (8) 

We seek a solution in expansion series as 



z = ^m'^^tp). 

p=0 



(9) 



a± 



w 
2" 



1± W1 + 



(10) 



One of the roots, Q3 , is unphysical, because it does not sat- 
isfy Eq. ^ at 0(m'^). By using all the 0th order roots in- 
cluding unphysical ones, /o is factorised as 

fo{z) = w*{w* - e*){z - ef{z - as){z - a+){z - a_). (11) 



4.2 1st order 

Next, we seek Ist-order roots. We put 2 = a± -I- mz^i-^ -\- 
Oim?). At the linear order in m, Eq. ((5]) becomes 



Z(i)fo{oL±) + fi{a±) = 0, 



(12) 



where the prime denotes the derivative with respect to 2. 
Thereby we obtain a Ist-order root as 



Ma±) 



(13) 



/o(a±) 

The similar manner cannot be applied to a case of e, because 
it is a doublet root with /o(e) = fo{^) ~ 0, while /q (e) 7^ 0. 
At 0(m2), Eq. (O can be factorised as 

(2(i)[(™*-e*)e+l]+e) 

X (2(i)[(e-io)(u>*e-H)-e]-fe(e-«;)) =0. (14) 
Hence, we obtain two roots as 



2(1) = 
2(1) = 



(e* — w*)e — 1 ' 
e(e — w) 



(15) 
(16) 



(e — w)(w*e + 1) — e 

Here, the latter root expressed by Eq. (|16p is unphysical and 
thus abandoned, because it doesn't satisfy the original lens 
equation ((2}. On the other hand, the former root by Eq. p5|l 
satisfies the equation and thus expresses a physically correct 
image. 

4.3 2nd Order 

First, we consider perturbations around zeroth-order solu- 
tions of a±. At 0{m?), Eq. (O is linear in 2(2) and thus 
easily solved for 2(2) as 



i(2) 



4)/o («±) + 22(i)/i(q±) + 2/2(a±) 
2/;(a±) 



(17) 



Next, we investigate a multiple root e. At O(m^), Eq. 
([5| becomes linear in 2(2). It is solved as 



2(2) = 



2fi)/o (6)+ 32^1) /i (e)+62(i)/2(6) 

6[2(i)/o(e) + /;(6)] 



(18) 



4.1 0th order 

At 0{mP), the lens equation becomes the fifth-order polyno- 
mial equation as /o = 0. Zeroth order solutions are obtained 
by solving this. All the solutions are e (doublet), 03 and a±, 
where we defined 



ct'i — 



4.4 3rd Order 



Around zeroth-order solutions of a±, Eq. ([Sjl at O(m^) is 



linear in 2(2) and thus solved as 



2(3) = 



/o("±) 



1 3 

[2(l)2:(2)/o (a±) + ^2(i)/o (a±) 



+2(2),/i(a±) + y(i)fiia±) + Hi)f2{a±)]- (19) 
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Also around the multiple root e, Eq. ((5]) at O(m^) be- 
comes linear in 2(2) . It is solved as 



1 



Z(3) 



^{i)/o'(f) + /;W 



24 



1 3 

+Z{i)Z(2)fi (e) + g2(i)/i (e) 
+2^(2) /2 (6)1, 



(20) 



where we used /2 (z) = 0. 

Table[T]shows a numerical example of perturbative roots 
and their convergence. 



5 PERTURBATIVE SOLUTIONS FOR A 
POLYNOMIAL FORMALISM 2: TRIPLET 
LENS 

In a binary case, we have only the single parameter m for 
the perturbations. For A'^ point masses, we have to take ac- 
count of couplings among several expansion parameters. In 
addition, the degree of the polynomial becomes A^^ + 1, so 
that we cannot write down the whole equation. In order to 
get hints for A'' point-mass lenses, in this section, we inves- 
tigate triple-mass lenses explicitly up to the second order in 
i>2 and Vs. 

The polynomial equation is rewritten as 



3 3 

J2 E(^2)"('^3)^-^<?fe)(p3)(^) = o, 



(21) 



where we defined 

S(0)(0)(2) = (2 - £2)^(2 - e3)^[(w* - £2)2 + 1] 

x[{w* — el)z + l]{w* — WW* z — w). (22) 

We seek a solution in expansion series as 



oc oc 



^=Yl Xl('^2)'"(j^3)'''2(p,)(p3). 



(23) 



P2=0P3=0 



5.1 0th order 



Zeroth order solutions are obtained by solving the tenth- 
order polynomial equation as g(o){o) — 0. The roots are £2 
(doublet), £3 (doublet), as, a4 and a±, where we defined 



Q3 



(24) 



For the same reason in the binary lens, 03 and 04 are un- 
physical, in the sense that it does not satisfy the lens equa- 
tion ([2|. By using all the 0th order roots, <?(o)(o) is factorised 
as 



g(0)(0)(z) = - £2)(lt)* £3)(2 - £2)^(2 - £3)^ 

x(2 — Cfi){z — ce4){z — a+)(2 — a_). 



(25) 



5.2 1st order 

Here, we seek Ist-order roots. The image position is ex- 
panded as 2 = Q± +/^22(i)(0) +!/32(0)(l) +0{ui, l'l,1^2V'i)- By 

making a replacement in notations as 2 <-» 3, one can con- 
struct 2(o){i) from 2(i)(o). Hence, we focus on 2(i)(o) below. 
At the linear order in V2, Eq. (O becomes 



2(i){0)ff(0)(0)(a±) + 3(i){o)(q±) = 0. 
Thereby we obtain a Ist-order root as 
5(i){0)(a±) 



Z(im - 



Horn 



(a±) 



(26) 



(27) 



For the triple mass lens system, the root £2 becomes 
triplet with g(oo)(e2) = 3(oo)('=2) = ff(oo)(<^2) = 0, while 
5(00) (^) 7^ 0- After rather lengthy but straightforward calcu- 
lations, Eq. ([5]) at 0(1^2) can be factorised as 

{Z(l)(0)[{w* ~ £2)£2 + 1] + £2) 
X {ziimiiw* ~ el)(.2 + 1] + £2) 

X (2(1){0)[(<;2 - w)(w*e2 + 1) - £2] + £2(£2 - w)) 

= 0. (28) 

Hence, we obtain three roots as 
£2 



Z(im - 



2{1){0) - 



Z{1){0) 



(£* -U)*)£2 - 1' 
£2 

{el -w*)e2 - 1' 

£2 (£2 — w) 



(29) 
(30) 

.^ -■ (31) 

(£2 — w){w''e2 + 1) — £2 

At the linear order in U2, true solutions for the triple lens 
system has to agree with that for the binary system, when 
one takes a limit as ^3 0. Therefore, out of the above 
three roots, ones expressed by Eqs. (|30p and (|3ip must be 
abandoned. 



6 PERTURBATIVE SOLUTIONS FOR A 
POLYNOMIAL FORMALISM 3: N 
POINT-MASS LENS 

In the previous section, we have learned couplings between 
the second and third masses. Now we are in a position to 
investigate a lens system consisting of N point masses. 
The polynomial lens equation ([5} is expanded as 



AT JV 



P2=0P3 = PN=0 



X5(P2){P3)--{Piv)(^) = 0- (32) 
For this equation, we seek a solution in expansion series as 



00 00 

EE-E 

P2=0P3=0 PN=0 
XZ(P2)(P3)---(Pn)- 



{^2^(1^37' ■ ■ ■ {i^nY 



(33) 



6.1 0th order 



Zeroth order solutions are obtained by solving the (A'^^ -f 
l)th-order polynomial equation as (/(o).- {o) ~ 0. The roots 



6 H. Asada 



Table 1. Example of perturbative roots in the single-complex-polynomial: We assume v = 0.1, e = 1 and two cases of ui = 2 (on-axis) 
and w = 1 + i (off-axis). In this example, the number of images is three. "Polynomial" in the table means all the five roots for the 
single-complex-polynomial. They are obtained by numerically solving the polynomial. Image positions are determined also by numerically 
solving the lens equation. They are listed in the column "Lens Eq.". In the table, "None" means that it does not exist. These tables show 
that, as we go to higher orders, the perturbative roots become closer to the correct ones for the single-complex-polynomial, including the 
two unphysical roots. 



Case 1 


On-axis 


u = 0.1 


e = 1 


w = 2 




Root 


1 


2 


3 


4 


5 


1st. 


2.43921 


-0.389214 


0.95 


-0.925 


0.975 


2nd. 


2.43855 


-0.388551 


0.95 


-0.924063 


0.974063 


3rd. 


2.43858 


-0.388519 


0.949938 


-0.924016 


0.974016 


Polynomial 


2.43858 


-0.388517 


0.949937 


-0.924013 


0.974013 


Lens Eq. 


2.43858 


-0.388517 


0.949937 


None 


None 



Case 2 


Off-axis 


u = 0.1 


e = 1 


w = 1 + i 




Root 


1 


2 


3 


4 


5 


1st. 
2nd. 

3rd. 


1.33716+1.40546 i 
1.33632+1.40363 i 
1.33634+1.40371 i 


-0.337158-0.355459 i 
-0.336316-0.354881 i 
-0.336275-0.354839 i 


0.95-0.05 i 
0.95-0.05 i 
0.95-0.05025 i 


0.025-0.925 i 
0.02625-0.9225 i 
0.0262813-0.922281 i 


0.975-0.025 i 
0.97375-0.02625 i 
0.973656-0.0263438 i 


Polynomial 


1.33633+1.40371 i 


-0.336272-0.354835 i 


0.950015-0.0502659 i 


0.0262762-0.922254 i 


0.973646-0.0263517 i 


Lens Eq. 


1.33633+1.40371 i 


-0.336272-0.354835 i 


0.950015-0.0502659 i 


None 


None 



((5]) at 0{v2) can be factorised as 

N 

H (^{i){0)-(0)[(^fe)'^2 + 1] + ea) 

k-2 

X (2;(i)(o)...{o)[™2(™*e2 + 1) + ^2] + e2W2) = 0. (37) 

By using this factorisation, we obtain TV roots. At the hnear 
order in 1/2, however, true solutions for the present lens sys- 
tem has to agree with that for the binary system, because 
one can take the limit as i/p — > for p > 3. Therefore, only 
the — e2/(u)2e2 + 1) out of the above A'^ roots is correct for 
the original lens equation. The same argument is true of any 



7 PERTURBATIVE SOLUTIONS FOR 

ZZ'-UV AL FORMALISM 1: BINARY LENS 

As shown above, an analytic polynomial formalism is appar- 
ently simple. When we solve perturbatively the polynomial 
equation, however, we find unphysical roots which satisfy the 
polynomial but does not the original lens equation. In the 
polynomial formalism, therefore, we are required to check 
every roots and then to pick up only the physical roots sat- 
isfying the original lens equation with discarding unphysical 
ones. It is even worse that the order of the polynomial grows 
rapidly as A'^'^ + 1, as the number of the lens objects increases. 
This means that the perturbative structure of the formalism 
becomes much more complicated as N increases. In this sec- 
tion, we thus investigate another formalism, which allows a 
more straightforward calculation especially without needing 
extra procedures such as deleting physically incorrect roots. 



are ai = —1/w*, a±, and ei (with multiplicity = A'^) for 
i — 2, ■ ■ ■ N , where for later convenience we denoted 

Wi =w - ei. (34) 

Like in the binary lens, m is unphysical, in the sense that it 
does not satisfy the lens equation ([2j. By using all the 0th 
order roots, (/(o) -.(o) is factorised as 

JV 

ff(0)...(0)(z) = {z ~ a+){z - a-)w*Y\_{i>'*) 

N N 

xl[(z-e,)''l[iz+±), (35) 



where this degree is + 1 in agreement with that of the 



polynomial equation. 



6.2 1st order 

Next, we seek Ist-order roots. In the similar manner in the 
double or triple mass case, we can obtain a Ist-order root as 

g(0)-(ifc)--(0)(«±) l'oR^ 
Z{0)..-(ifc)--{0) -, 7—; — , (,3b) 

where Ik denotes that the fc-th index is the unity, namely 
Pk = 1- 

For A'^ point mass lens systems, a root efe is multiplet 
with multiplicity = N. Without loss of generality, we choose 
£2 as a root in the following discussion. Calculations done 
above for a double or triple mass system suggest that Eq. 
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First, we focus on a binary case for its simplicity. The 
lens equation is rewritten as 



C{z, z*) = mD{z*), 
where we defined 



C(z,z*) = w 
D{z*) 



■z + - 



1 1 

z* z* — 



(38) 

(39) 
(40) 



One of advantages in this zz'-formulation is that the 
master equation (|38p is linear in m. Therefore, counting or- 
ders in m can be drastically simplified when we perform 
iterative calculations. On the other hand, an analytic poly- 
nomial is second order in m. In fact, practical perturbative 
computations in the polynomial formalism are quite compli- 
cated, in the sense that several different terms (/o, /i and /2 
for a binary case) may make the same order-of-magnitude 
contributions at each iteration step. 

We seek a solution in expansion series as 



z = y m z^k)- 

k=0 



The complex conjugate of this becomes 

oo 

z* = y^m'^z'fc). 



(41) 



(42) 



According to these power-series expansions of z and z* , both 
sides of the lens equation are expanded as 



oc 



(43) 



(44) 



where Cj^) and D(t,) are independent of m. At 0{m''), Eq. 
(1381) becomes 



(45) 



C(fc) — D(^k_-i), 



which shows clearly a much simpler structure than a polyno- 
mial case such as Eqs. I|17|l and (|18|l . Equation (|4U|I indicates 
that D{z*) has a pole at z* = e* . Therefore, we shall discuss 
two cases of 2(o) 7^ e and Z(o) = e, separately. 



The L. H. S. is purely real so that the R. H. S. must be 
real. Unless w = Q, therefore, one can put 2(o) = by 
introducing a certain real number A. By substituting 2(o) = 
Aw into Eq. (|47|) . one obtains a quadratic equation for A as 



WW* A — WW* A - 
This is solved as 



1 = 0. 



A 



1 



1± W1 + 



(48) 



(49) 



which gives the Oth-order solution. 

In the special case of w = 0, Eq. (I47|l becomes |2(o)| = 
1, which is the Einstein ring. In the following, we assume 
w^G. 

7.2 1st order (^(o) 7^ e) 

In units of ^(o) , the expansion series of z is normalised as 



z = 2(0) y^m'°a(fc), 



(50) 



where we defined CT(fc) = z^^^jz^Qy 

First, we investigate a case of z(o) 7^ e. At the linear 
order in m, Eq. psp becomes 



2(1) + 



'(1) 



+ 



(51) 



'(0)'' "^(0) -^(0) 

In order to solve Eq. (|51[) . we consider an equation linear 
in both z and z* as 



(52) 



(53) 



2 + az* = 6, 

for two complex constants a,b £ C. 

Unless \a\ = 1, the general root for this equation is 

b — ah* 

z = . 

1 — aa* 

This can be verified by a direct substitution of Eq. ((53} into 
Eq. (f52)) . If \a\ = 1, Eq. ((52} is underdetermined, in the 
sense that it could not provide the unique root without any 
additional constraint condition on 2 and 2*. 

By using directly Eq. (|53|) . Eq. (|51|l is solved as 



1 



^{0)(^(0))^ ^ 



* 2 * 
^ ^(0)^(0) 



(0) 



HO) 



HO) 



(54) 



T.l 0th order (2(0) 7^ e) 

Zeroth order solutions are obtained by solving the equation 
as 



C(2{0), 



= 0. 



(46) 



The solution for this is the well-known roots for a single 
mass lens. In order to help readers to understand the 22*- 
dual formulation, we shall derive the roots by keeping both 
2 and 2*. In conventional treatments, a single lens case is 
reduced to one-dimensional one by choosing the source di- 
rection along the x-axis in vector formulations or the real 
axis in complex ones. Eq. (|46p is rewritten as 



HO)- 



(47) 



7.3 2nd order (2(0) / e) 

At 0(m2), Eq. dSHI is 

2(2) + a2 2(2) = 62, 

where we defined 
1 



Cl2 



62 



(-(0))^' 



-Ad + 



«(0) 



Here, D(i) is written as 



(1) 



^(*0) 



+ 



'(1) 



(55) 

(56) 
(57) 

(58) 



'(0) 



8 H. Asada 



Z(2) = 



By using the relation (|53p . Eq. (|55|) is solved as 

62 " £12 &2 



(59) 



7.4 3rd order and nth order (2(0) 7^ e) 

Computations at 0{m?) are similar to those at 0{m?) as 
shown below. At O(m^), Eq. (|38p takes a form as 



Z{3) + as 2*3) = 63, 

where we defined 
f 



013 = 



63 = 



I>(2)+ 



'(0) 



Here, _D(2) is written as 



D 



(2) 



2(3) 



'{2) 



('^(1)) 



+ 



<5(2) 



(0) 



^(0) (^(0) 

Using the relation (|53|l for Eq. (|60|l . we obtain 

_ 63 - asbs 
1 — aza% 



(60) 

(61) 
(62) 

(63) 
(64) 



In the similar manner, one can obtain iteratively nth- 
order roots 2(„) , which obeys an equation in the form of Eq. 
(|52p . and thus can use Eq. ((53} to obtain 



7.5 0th and 1st order (2 



(0) 



Next, we investigate the vicinity of z = e, which is a pole of 
D. The other pole of D is 2 = 0, which makes also C{z, 2*) 
divergent. Therefore, 2 = and its neighbourhood are aban- 
doned. Let us focus on a root around 2 = e. 

We assume z — e-\- mz(^i-^ + 0{m?). Then, the relevant 
terms in expansion series of C and D become 



(0) 



D 



(-1) 



'(1) 



(65) 
(66) 



where the index —1 means that the inverse of m appears be- 
cause of the pole at e. Therefore, the lens equation at 0{mP) 
becomes linear in 2(\j without including 2(i). Immediately, 
it determines 2(\j. Its complex conjugate becomes 



2(1) = 



(w* — e*)e + 1 



(67) 



This shows a clear difference between 2(o) = e and 2(o) 7^ e 
cases. Equation (|51|l for the latter case contains both 2(i) 
and 2(*i), so that we must use a relation such as Eq. (f53)) . 



7.6 2nd, 3rd and nth order (2 



(0) 



Next, we consider the lens equation at O(m^), namely 

(68) 



C(i) = D(^()y This determines 2^2) ^ 



<2) = (2(1))' (C(i) - ^) 
where we may use 







-1 







1 

X 



Figure 2. Pcrturbative image positions for a binary lens case. 
This plot corresponds to Tables [T] and |2] The lenses (ei = 0, e2 = 
1) and sources (ui = 2 and ui = are denoted by filled squares. 

The image positions are denoted by filled disks. Pcrturbative im- 
ages at the 1st, 2nd and 3rd orders are overlapped so that we 
cannot distinguish them in this figure. 



(1) 



-2(1) 



^(1) 
(,,)2 



(69) 



Let us consider O(m^) to look for 2(3). Equation of 



(2) 



provides 2^3^ as 



, _ (* n2^ 
2(3) - (,2(1) j 0(2) + -J^Z^ 



+ 



"(1) 



where we can use 



C(2) — —2(2) — 



(2) 



'(2) 



(2*1 



(70) 



(71) 



By the same way, one can obtain perturbatively nth-order 
solutions 2(„) around 2(oi = e. 



(0) 



Table [2] shows an example of perturbative roots in 
the dual-complex- variables formalism and their convergence. 
Tables [T] and [2] suggest that the polynomial approach and the 
dual-complex-variables formalism are consistent with each 
other, regarding the true images. Figure [2] shows image po- 
sitions on the lens plane, corresponding to these tables. 



8 PERTURBATIVE SOLUTIONS FOR 

ZZ*-DUAL FORMALISM 2: LENSING BY N 
POINT MASS 

The purpose of this section is to extend the proposed method 
to a general case of gravitational lensing by an arbitrary 
number of point masses. 

The lens equation is written as 



C(2, 



JV 

E 



VkDk{z*), 



(72) 
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Table 2. Example of perturbative images via the dual-complex-variables formalism: We assume the same values for parameters as Table[T] 
Good agreements with these tables suggest a consistency between the single-complex-polynomial and the dual-complex- variables formalism, 
regarding the true images except for unphysical roots. 



Case 1 (On-axis) 


u = 0.1 


e = 1 


w = 2 


Root 


1 


2 


3 


1st. 


2.43921 


-0.389214 


0.95 


2nd. 


2.43855 


-0.388551 


0.95 


3rd. 


2.43858 


-0.388519 


0.949938 


Lens Eq. 


2.43858 


-0.388517 


0.949937 



Case 2 (Off-axis) 


u = 0.1 


e = 1 


w = 1 + i 


Root 


1 


2 


3 


1st. 

2nd. 

3rd. 


1.33716-1-1.40546 i 
1.33632-1-1.40363 i 
1.33634-1-1.40371 i 


-0.337158-0.355459 i 
-0.336316-0.354881 i 
-0.336275-0.354839 i 


0.95-0.05 i 
0.95-0.05 i 
0.95-0.05025 i 


Lens Eq. 


1.33633-1-1.40371 i 


-0.336272-0.354835 i 


0.950015-0.0502659 i 



where C{z, z*) was defined by Eq. (|39|) and we defined 



z' z' ~ el 



(73) 



C{z,z*) and Dk{z*) in the lens equation (|72p are ex- 
panded as 



oc oc oo 



c{z, z*) = E E • ■ ■ E (^^2)'' ^"'^"^ ■ ■ ■ ^"^y 

P2=0P3=0 PJV=0 
^C'(P2)(P3)-"{Pn)(^' ^ )' 



CXD OO OO 



P2=0P3 = PJV=0 



where C, 



fe(P2)(P3)---(Pjv) 

and Di- 



(74) 



(75) 



^) are independent 



(P2){P3)'--(P«) "^na ^k(p2)(p3)--{PN 

of any Vi. We seek a solution in expansion series as 

oo oo oo 
P2=0P3=0 PJV=0 



where z, 



(P2)(P3)---{Pn) 



is a constant to be determined itera- 



tively. The perturbed roots are normalised by the zeroth- 
order one as 



^(P2)(P3)'--(PW) 
2(0)--(0) 



(77) 



Equation (|73p shows that Dk{z*) has a pole at z* = e^. 
Therefore, we shall discuss two cases of ^(o) 7^ or 2(o) — et, 
separately. 



8.1 0th order (z(o) (o) 7^ for * = li ' ■ ■ i ^) 

Zeroth order solutions are obtained by solving the equation 
as 



C(2,2*) = 0. 



(78) 



This has been solved for the binary lens case. The solution 
is given as 



= Aw, 



with the coefficient A defined by Eq. (|49 



(79) 



8.2 1st order (z(o)...(o) 7^ U for i = 1, ■ ■ ■ , A*') 
At the linear order in Vk, Eq. (|72|l is 

C(o)...{ij.)...(o) = J^fc-Dfe(o)...{o), (80) 

where 1^ denotes that the fc-th index is the unity. This equa- 
tion is rewritten as 

Z{0)..-(ifc).-{0)+a{0)-.{ifc)...{0) x4o)-.(it3)...{o) = ''{o)...(ifc)--(o), (81) 
where we defined 

1 



a(o).-{ifc)...(o) 

&(0)--(lfc)--(0) 



(^('o)--{o))^' 



(82) 
(83) 



By using Eq. (|53|l . we obtain 



2(0)...(lfc)...(0) 

&(0)..-(lfc)...(0) - "{0)..-(lfc)..-(0)6{0)...(lfc)...(0) 



1 - a 



(()).. .(lfc)...(0)ffl(o)...(ij^)...(0) 



(84) 



8.3 2nd order (z(o)...(o) 7^ for i = 1, ■ ■ • , Af) 

Let us consider two types of second-order solutions as 
2(0)---(2fc)...(o) and 2(o)...(i^)...(ij)...(o) for k ^ £, separately. 

First, we shall seek 2:(o)...(2fc)...(o) • At 0{iyl), Eq. ((72} 
becomes 



2(0)-.(2fe)...(0) + a(0)...(2fc)...(0)2:(0)-.(2fc)-.(0) 
= &(0)--(2fc)-.(0), 



where we defined 
a(o)--(2fc)-.(o) = 



(^('o).-{o))^ 



(85) 



(86) 
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3(0)--{2fc)...(0) 



= -D 



where 



(0).-(lfc)...{0) 



fc(0)---(lfc)---(0) 



is written as 



+ 



( °"(*o)--(ifc)--(o))^ 

^(*0)--(0) 



(87) 



D 



fc(0)...(lfc)...(0) 



2(0)--(lfc)--(0) 



(^(*0)...{0))^ 



(O)--(O) 



-tfc 



By using the relation (I53|l for Eq. (|85p . we obtain 



^(0)--(2fc).-{0) 

&(0)...{2fc)...(0) - "(0)...(2;,)...(0)6{0)-.(2t_)'--(0) 



1 - a 



(0). ..(2fc)... (0)^(0). ..(2fc)...(0) 



(89) 



Next, let us determine Z(o)...(ij.)...(if )...(o) • At 0{vkvt) 



(90) 



for k <£, Eq. ((72} becomes 

^(0)--(lfc)--(lf)--(0) 

+"(0)--(lfe)--(l«)--(0)Z(*0)--(lfc)--(lf)--(0) 



= 6 



(0)--(lfc)--(lf)--(0), 



where we defined 

a(0). ..(ifc).. ■(!«). ..(0) 



(^{*o).--(0))^ 



(91) 



6(o)...(ifc)..-(if)...(o) - --Dfe(o)...(if)...(o) - -Df(o)...(i^)...(o) 



'-(O)--(O) 



Here, -Dfc(o)-.-(if )...(o) and (o)...(i^)...(o) 



are written as 



-Dfc(0).-(if).-{0) 



2(0)-.(lf)-.(0) 



{O)-.(O) 



f(0)--(lfc)--(0) 



-Z(0)--(lfc)--(0) 



(^{*0)---(0))^ (^('o)...(0) '''^)^ 



,(93) 



.(94) 



By using the relation (|53p for Eq. (|9Up , we obtain 

^(0)--(lfe)--(lf)--(0) 

_ &(0)...(lfe)...(lf)...(0) - «(0)...(lfc)...(lf)-.(0)&*0)-.(lfc)-.(lf)...(0) 



1 - a 



{0)..-(lfc).. -(If).. -(0)0(0). ..(ij^)...(lj)...(0) 



7? 



(0)--{-lfc)--{Q) 



^(0)--(lfc)--{0) 



(98) 



where —Is, means that the inverse of i^k appears because of 
the pole at efe. Therefore, the lens equation at 0{vk) becomes 
linear in z*o)...(i^)...(o) without including 2;(o)...(ij^)...(o) . Im- 
mediately, it determines 2:*q)...(i^)...(o)- Hence, its complex 
conjugate provides 

(it)* - e^)efc + 1 
8.5 2nd order (z(o)...(o) = tfc) 

Here, we consider the lens equation at 0{iyl), namely 
C(o)...{ifc)...(o) = -D(o)...(o), where 



■^(0)--{0) - -T + 73 



(0)--{2fc)...(0) 



(100) 



{0)--(lfc).--{0)^ 



Hence, we obtain z, 
jugate as 



(0).--(2^,)-.{0) 



and thereby its complex con- 



Z{0).--(2fc).-{0) - l2(0).-{lfc)...(0)J 



where C(o)...(i^)...(o) becomes 



(0)--(lfc)--(0) 



(0)--(lfc)--(0) 



2(0)--(lfc)--(0) + 



^(0)--(lfc)--(0) 



77 , (101) 



(102) 



Next, we consider a root at 0(v\vl), where we can as- 
sume k < £ without loss of generality. At this order, the 
inverse of fk appears. The lens equation at 0{u}) becomes 

C(o)...{ij)...{o) = -Dfc{o)...(-i^.)...(o) + -D«{o)...(o), (103) 



where 



D 



fe(o)---(-ifc)---(if)---(o) - (• 
Hence, we obtain z 



^(0)--(u).--(i^).--(0) 



(104) 



plex conjugate as 

Z{0)...(lfc)...{l^)...{0) 
= (2(0)--(lfc).--(0))^ 



(0)--{lfc)--{0)^ 

(o)...(ifc)...{i,)...o) and thereby its com- 



(95) 



{CI, 



(o)--(if)--(o) - -D£(o)...(o)j , 



8.4 0th and 1st order (2(o)...(o) = efc) 

Next, we investigate the vicinity of z — e^, which is a pole 
of Dk- The other pole of Dt is ^ = 0, which makes C{z, z*) 
divergent. Therefore, z — and its neighbourhood are aban- 
doned. Let us focus on a root around 

2(0)---(0) = Efc- (96) 

If we admitted .Z(o)...(if ) -.(o) around for I k, only 
the Dk function would contain the inverse of vi, which 
introduces a term at Oivk/vi) in the lens equation and 
leads to inconsistency. Namely, the lens equation prohibits 
2(o) - (i(.) - (o) around tfe for I / k. This agrees with the poly- 
nomial case. We thus assume z = et + i/fc2:(o)...(ij.)...(o) + 
0{ul). Then, we obtain 



where C(o)...{i,)...(o) and D, 



are written as 



(0)--(lf)--(0) 



-2:(0)...(i^)...(0) 
1 1 



e{0)-{0) 



(105) 



, (106) 
(107) 



C(o)...(o) = w-ek + —, 



(97) 



This direct computation shows that Z(o)...(ij,)...{if )...{o) does 
not exist because Z{o). ■(i;)...(o) is prohibited in the vicinity 
of ek ■ This is also consistent with the polynomial case at the 
second order. 



8.6 Magnifications 

Before closing this paper, it is worthwhile to mention mag- 
nifications by A'^ point-mass lensing in the framework of the 
present perturbation theory that is intended to solve the 
lens equation to obtain image positions. The amplification 



Gravitational Lensing by N Point Mass 11 



< 




< 




-0.5 0.0 0.5 
t 



1.0 



0.0004 



0.0002 



< 0.0000 



-0.0002 



-0.0004 




-2-10 1 
t 



0.15 



0.10 



^ 0.05 



0.00 




Figure 3. Light curves by two methods. One is based on a 
numerical case that the lens equation is solved numerically. The 
other is due to the first order approximation. The top figure shows 
that the two curves are overlapped, where A denotes the total 
amplification. The bottom panel shows the residual by the two 
methods. The residual is defined as the difference between A com- 
puted numerically and A in the linear approximation. We assume 
the source trajectory a.sw = lA + it. Here, the time t is in units of 
the Einstein cross time, which is defined as 8e/v±^ for the trans- 
verse angular relative velocity. The lens parameters aie 1^2 =0.1 
and e = 1. 



factor is the inverse of the Jacobian for the lens mapping. It 
is expressed as 



\de) 



Figure 4. Light curves by two methods. In this figure, we assume 
a different source trajectory as ui = 0.8 + it. The lens parameters 
are the same as V2 =0.1 and e = 1 in Fig. |3] The solid curve in 
the top panel denotes a case when the lens equation is numeri- 
cally solved. The dotted curve is drawn by using the linear order 
approximation. The bottom panel shows the residual between the 
two curves. 



d{w, 


w*) 


-1 




d{z, 




) 




dw 


2 


dw 


r 


dz 




dz* 





(108) 



where the terms in the last line can be computed directly 
by a derivative of Eq. Q, the lens equation in a complex 
notation. Amplifications of each image are obtained by sub- 
stituting its image position into Eq. (|108|l . Practical numer- 
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dz_ ^ 9z(o)...(o) ^ ^,^ ^^(o)--(ifc)--(o) 
9-1(7 dw ^—^ dw 

k 

and for self-interacting images 



dz 
dw 



dz 



Vk- 



(0) --(lfc)--(0) 



(109) 



(110) 



where we used that is a constant. 

For the simplicity, we assume Vk = 0{1/N) for a large 
case. Then, the linear order term in self-interacting images is 
0{1/N), and thus they become negligible as A'' ^ oo. On the 
other hand, mutually-interacting ones have non-vanishing 
terms even at the zeroth order. Hence, they can play a cru- 
cial role. 

However, we should take account of a spatial distribu- 
tion of lens objects. If they are clustering and thus dense 
at a certain region, then the total flux of light through 
such a dense region is not negligible any more. Let us de- 
note the fraction of the clustering particles by /. Total 
contributions from such clustering self-interacting images 
are estimated approximately as a typical image magnifi- 
cation multiplied by the number of the particles, namely 
fN X !^fc(~ 1/-^) ~ 0{f), which does not vanish even as 
N ^ oo. Figures |6] and [7] show an example of a large A'^ 
case, where A'' is chosen as 1000. 



Figure 5. Graph representations of interactions among point 
masses for images at the second order level. The top and bot- 
tom graphs represent a mutually-interacting image and a self- 
interacting one, respectively. 



ical estimations may follow this procedure. For illustrating 
this. Figs. [3] and U show examples of light curves by a binary 
lens via the perturbative approach. These curves are well 
reproduced. However, double peaks due to caustic crossings 
cannot be reproduced by the present method. 

As an approach enabling a simpler argument before go- 
ing to numerical estimations, we use the functional form 
of perturbed image positions. In the perturbation theory, 
lensed images can be split into two groups. One is that 
their zeroth-order root is not located at a lens object 
(2(0)... (0) 7^ £*:)• In the other group, zeroth-order roots orig- 
inate from a lens position at et- We call the former and 
latter ones mutually-interacting and self -interacting images, 
respectively, because all the lens objects make contributions 
to mutually-interacting images at the linear order as shown 
by Eq. (|84p . On the other hand, self-interacting images are 
influenced only by the nearest lens object at at the lin- 
ear and even at the second orders as shown by Eqs. I|99p 
and (|101l) . Figure [5] shows graph representations for the two 
groups of images. 

For the simplicity, we consider stretching of images 
roughly as \dz/dw\, though rigorously speaking it must be 
the amplification. Table 1 and Equation (|76|) mean that the 
complex derivative becomes for mutually-interacting images 



9 CONCLUSION 

Under a small mass-ratio approximation, this paper devel- 
oped a perturbation theory of N coplanar (in the thin lens 
approximation) point-mass gravitational lens systems with- 
out symmetries on a plane. The system can be separated 
into a single mass lens as a background and its perturbation 
due to the remaining point masses. 

First, we investigated perturbative structures of the 
single-complex-variable polynomial, into which the lens 
equation is embedded. Some of zeroth-order roots of the 
polynomial do not satisfy the lens equation and thus are 
unphysical. This appearance of correct but unphysical roots 
is consistent with the earlier work on a theorem on the max- 
imum number of lensed images (Rhie 2001, 2003). However, 
the theorem never tells which roots are physical (or unphys- 
ical). What we did is that unphysical roots are identified. 

Next, we re-examined the lens equation in the dual- 
complex-variables formalism to avoid inclusions of unphysi- 
cal roots. We presented an explicit form of perturbed image 
positions as a function of source and lens positions. As a 
key tool for perturbative computations, Eq. (|53|) was also 
found. For readers' convenience, the perturbative roots are 
listed in Table [3] If one wishes to go to higher orders, our 
method will enable one to easily use computer algebra soft- 
wares such as MAPLE and MATHEMATICA. This is be- 
cause it requires simpler algebra (only the four basic opera- 
tions of arithmetic), compared with vector forms which need 
extra operations such as inner and outer products. 

There are numerous possible applications along the 
course of the perturbation theory of N point-mass gravi- 
tational lens systems. For instance, it will be interesting to 
study lensing properties such as magnifications by using the 
functional form of image positions. Furthermore, the validity 
of the present result may be limited in the weak fleld regions. 



Gravitational Lensing by N Point Mass 13 



1.0 
0.5 
0.0 
-0.5 
-1.0 




-1.0 -0.5 0.0 0.5 1.0 

X 



1.0 
0.5 
0.0 
-0.5 
-1.0 



-1.0 -0.5 



0.0 

X 



0.5 1.0 



1.0 
0.5 
0.0 
-0.5 
-1.0 



• . : .•• •.••••» •.i'". ». . . 



-1.0 -0.5 0.0 0.5 1.0 

X 



Figure 6. Example of a large N case. Here, we assume a trun- 
cated isothermal sphere projected onto a single lens plane with 
A'^ = 1000, where the truncation radius is the unity. For the sim- 
plicity, we assume equal masses. The source located at 0.25 is 
denoted by the circle. The top figure shows locations of the N 
point masses on the lens plane. The bottom shows a plot of im- 
age positions by using the perturbative solutions at the second 
order. In practice, the linear-order and second-order roots make 
no difference distinguishable by eyes in the figure. 



It is important also to extend the perturbation theory to a 
domain near the strong field. 

Our method considers only the images which exist in 
the small mass limit as i^i 0. The number of the im- 
ages that admit the small mass-ratio limit is less than the 
maximum number. This suggests that the other images do 
not have the small mass limit. Therefore, it is conjectured 
that positions of the extra images could not be expressed 



Figure 7. Plot of image positions with lensing amplification 
for a case of = 1000. The source and lenses are the same as 
those in Fig. [6] Here, we take account of amplifications by lensing. 
The area of a disk corresponding to each image is proportional to 
the magnification factor in arbitrary units. Large amplifications 
near ±1 are caused by the mutually-interacting images. On the 
other hand, a concentration of small but many images around 
the center are due to the self-interacting images, because lens 
objects have a large number density there. These three regions 
may correspond to three images for a singular isosphere lens in 
the limit of Af ^ oo. 
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Figure 8. Einstein ring broken by the lens discreteness due to the 
finite- effect. The lenses are the same as those in Figs. [6] and [7] 
The source is located at the origin of the coordinates. Amplifica- 
tions are taken into account. The area denotes the magnification 
factor in arbitrary units. 
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Table 3. List of the coefficients in perturbative positions of im- 
ages lensed by N point masses: The image positions are expressed 
in the form of 2: = ' ' ' Ep„ (^2)"' ' ' ' ('^iv)''"^{P2)--(Piv)- 

The top and bottom panels show the cases of 2(o) - (0) = ^rid 
2(o)---(o) 7^ ^si respectively. In the columns, "None" means that 
the corresponding coefficient does not exist. 

Case 1: ^ (i = 1, ■ • ■ , A^) 

2(0). ..(0) Eq. dZHJ 

2(o)...(i^)...(0) Eqs. (ESJ-lHl 

2(o)...(2^)...(0) Eqs. (|86]l-l[89ll 

2(o)...{i^)...(ij)...(0) Eqs. (|9lJ-ll95ll 

Case 2: ^(O)---(O) = ^fc 

2(0)... Eq. (isg 

2{0)---(lfc).--(0) Eqs. (I99J 

2(0)---(lf)---(0) for e^k None 

H0)-(2^)--(0) Eq. JlOlll 

^(0)---(lfc)---(if)---(0) for £ 7^ A; None 

as Maclaurin series in mass ratios. This may be implied also 
by previous works. For instance, the appearance of the max- 
imum number of images for a binary lens requires a finite 
mass ratio and the caustic crossing (Schneider and Weiss 
1986). Regarding this point, further studies will be needed 
to determine positions of all the images with the maximum 
number as a function of lens and source parameters. 
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